Could Z c (3900) be a I G J P = 1+1+ D*D molecular state? 
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^, ', The recently observed narrow resonance Z c (3900) is examined with the assump- 

tions as a D*D molecular state with quantum numbers I J = 1 + 1 + . Using QCD 
sum rules (QCDSR), we consider contributions up to dimension six in the opera- 
tor product expansion (OPE) at the leading order in a s . The mass is obtained as 

<C ; (3.91 ± 0.19) GeV which coincide with the Z c (3900). 

PACS numbers: ll.55.Hx, 12.38.Lg, 12.39.Mk 
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The charmonium-like exotic states are of great interest as they provide a satisfied window 
both in studying the dynamics of light quarks interact with heavy quarks and testing the 
standard model itself. In the past years, so many exotic states have been observed by the 
collaborations such as CDF , DO, BELLE and BABAR. They are also investigated with 
various models theoretically. These exotic states can be illustrated either as molecular 
states or multi-quark states or hybrids(see reviews ll|-p] and references therein). Although 
there is even no one conformed conclusion on these configurations, it is still interesting to 
investigate the inner structure of these states. 

Recently, the BEC Collaboration [7| reported a new enhancement structure in the 
Tt^J/ty invariant mass spectrum of the Y(4260) — > J/^tt + it~ decay. The mass and width 
of this state is M = (3899.0 ± 3.6 ± 4.9) MeV/c 2 and T = (46 ± 10 ± 20) MeV/c 2 . 
BELLE confirm this observation with mass M = (3894.5 ± 6.6 ± 4.5) MeV/c 2 and width 
T = (63 ±24 ±26) MeV/c 2 , where the errors are statistical and systematic, respectively 8j. 
Before this observation, Its existence is predicted by the ISPE mechanism [9| and in the 
tetraquark scheme [lOJ]. Many attempts have been made to investigate its possible con- 
figuration with various models, most of which support the D*D molecular structure with 
J p = 1+ In Ref. [ll|, They give an explanation of Y(4260) as DD 1 (2420) + DD t (2A20) 



X 



and inteprete Z c (3900) as D*D bound state. Based on heavy quark spin symmetry and 
heavy flavour symmetry, Guo predict Z c (3900) as the isovector D*D partners of the 
Zb(10610) [12|. In Ref. 14j, the authors discuss the Z c (3900) with quantum number 



J PG = 1 ++ within tetraquark model and molecular model, and investigate its decay modes. 



In Ref. 15], Z c (3900) is studied in the molecular model as well as the hadro-charmonium 



and tetraquark schemes, and by the ways of distinguishing by further experimental studies. 



In Ref. [16], it is argued that the molecular interpretation is less likely and this could be 
investigated with QCD sum rules. 

Due to the asymptotic property of the QCD, study of the hadron physics have to concern 
about the nonperturbative effect which is difficult in quantum field theory. There are many 



method to estimate the mass of a hadron, among which QCD sum rule(QCDSR) 17H21J] is a 
fairly reliable one. Quantum numbers compatible with the experiment are the fundamental 
ingredients in QCDSR analysis of composite particles. Since Z c (3900) was observed in the 
F(4260) — > Z c (3900)7r decay process, the assignment of I G = 1 + is known. In this article, 
by assuming Z c (3900) as a D~D*° molecular state with I G J P = 1 + 1 + , we investigate the 
mass of this possible molecular configuration within the framework of QCD sum rules. We 
construct the following current to interpolate the Z c (3900) molecular state 

f = -^[(m7 5 c)(c7^) + (u^c)(ci-f 5 d)}. (1) 

The rest of the paper is organized as three parts. The QCDSR for the Z c (3900) is 
derived in Sec. [Ill with contributions up to dimension six in the OPE. The numerical 
analysis is presented to extract the hadronic mass at the end of this section. Sec. IIHI is 
the summary and conclusion. 

II. QCD SUM RULES FOR Z c (3900) 

In the QCDSR approach, the mass of the particle can be determined by considering the 
two-point correlation function 

W{q 2 ) = i [ d 4 xe iq - x (0\T[f{x)f + {0)}\0). (2) 



Lorentz covariance implies that the two-point correlation function can be generally param- 
eterized as 

it-V) = (^r - <T)n (1) (g 2 ) + ^n(°V). (3) 

gZ gZ 

The term proportional to g^ v will be chosen to extract the mass sum rule, since it gets 
contributions only from the 1 + state. The QCD sum rule attempts to link the hadron 
phenomenology with the interactions of quarks and gluons. It contains three main ingredi- 
ents: an approximate description of the correlation function in terms of intermediate states 
through the dispersion relation, a description of the same correlation function in terms of 
QCD degrees of freedom via an OPE, and a procedure for matching these two descriptions 
and extracting the parameters that characterize the hadronic state of interest. 

In the phenomenological side, the correlation function is calculated by inserting a com- 
plete set of intermediate states. Parameterizing the coupling of the Z state to the current 



j M as 



(0\j,\Z) = Ae M . (4) 



U^(q 2 ) can be expressed as 



n() (g 2 ) = T7^-3 + r/ rfs 1 ^ . ( 5 ) 
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where Mz denotes the mass of the molecular state, and so is the threshold parameter. 
In the OPE side, UP^(q 2 ) can be written as 

IT«(g 2 )=/ &£— ^, (6) 

J4m2 S — g 2 

where the spectral density is p OPE (s) = -Imllv i /(s). Making quark-hadron duality as- 
sumption and a Borel transformation, the sum rule is obtained by matching the two sides: 

X 2 e -Ml/M* = f S ° dsp OPE {s)e -s/M^ (7) 

J4m| 

with M 2 the Borel parameter. 

In the OPE side, we work at the leading order in a s and consider vacuum condensates 



up to dimension six, with the similar techniques in Refs. 22|]. In order to consider the 
isospin violation, we keep the terms which are linear in the light-quark masses m u and 
nid- After some tedious OPE calculations, the concrete forms of spectral densities can be 
derived: 



p OPE (s) = pP ert (s) + p {qq) (s) + ptf & \s) + p {9 ^ Gq \s) + p^\s) + p {gZG3) (s), (8) 



with 
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Pert (s) = 3 r- * f-° | (1 _ a _ p)(1 + a + m , nc] s y 
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iaa\, s 3 (gg) ra max da p- a d/3 . _., n _. . , 2 

P W (s) = -^Tmc -j/ -^(« + /3)(l + « + /3)r(m c ,s) 2 
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(g 3 G 3 ) 2 r«m« /•!-« rf/3 



p(s . G , (3) = |o) m , r— da r- <v (1 _ Q _ W(1 + Q + w 



B B 3 
(g 3 G 3 ) r^rnax fl-a rf/3 _ 

013^6 / j R ~R3 ' 

A /I J Oimin J Pmin t J 



+ ini da JL(i- a -/3)(l + a + f3)r(m c , s), (9) 



with r(m c , s) = (a + /3)m 2 — a/3s. The integration limits are given by a m i n = (1 — 

yl -4m2/ s )/2, a max = (l + \Jl- 4m%/s)/2, and /3 min = am;?/(sa - m 2 c ). 

Taking the derivative of Eq. (|7|) with respect to -jk^ and then dividing by itself, we arrive 
at the mass of the molecular state 

M 2 Z = r dsp OPE se- s l M2 / r ds P OPE e- s ' M2 . (10) 

Before the numerical analysis of the equation f jTOj) . we first specify the input parameters. 
The quark masses are taken as m u = 2.3 MeV, uid = 6.4 MeV, and m c = (1.23 ±0.06) GeV 
23| . The condensates are (uu) = (dd) = (qq) = -(0.23±0.03) 3 GeV 3 , (gqa-Gq) = m 2 (qq), 
m 2 = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 , and (g 3 G 3 ) = 0.045 GeV 6 0. Complying with the 
standard procedure of the sum rule, the threshold Sq and Borel parameter M 2 are varied to 
find the optimal stability window. There are two criteria (pole dominance and convergence 
of the OPE) for choosing the Borel parameter M 2 and threshold Sq. 

The contributions from the high dimension vacuum condensates in the OPE are shown 
in FigJTJ We have used ^/sq > 4.2 GeV. From this figure it can be seen that for M 2 > 
2.50 GeV 2 , the contribution of the dimension-6 condensate is less than 20% of the total 
contribution and the contribution of the dimension-5 condensate is less than 26% of the 
total contribution, which indicate a good Borel convergence. Therefore, we fix the uniform 
lower value of M 2 in the sum rule window as M^ in = 2.50 GeV 2 . The upper limit of M 2 
is determined by imposing that the pole contribution should be larger than continuum 
contribution. FigfJ] shows that the contributions from the pole terms with variation of the 
Borel parameter M 2 . We show in Table [I] the values of M^ ax for several values of v/so. 
In Figjnj we show the molecular state mass, for different values of y/so, in the relevant 
sum rule window. It can be seen that the mass is stable in the Borel window with the 
corresponding threshold y/so- The final estimate of the I G J P = 1 + 1 + molecular state is 



TABLE I: Upper limits in the Borel window for the I J p = 1 + 1 + D*D current obtained from 
the sum rule for different values of y/so- 



V^ (GeV) 


Mi ax (GeY 2 ) 


4.2 


2.8 


4.3 


2.9 


4.4 
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4.5 
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3.4 
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FIG. 1: The OPE convergence for the molecular state. The contributions from different terms 
with variation of the Borel parameter M 2 in the OPE. The a and b correspond to the contributions 
from the D = 6 term and the D = 5 term, respectively. The notations a, /3, 7, A and p correspond 
to the threshold parameters ^/sq = 4.2 GeV, 4.3 GeV, 4.4 GeV, 4.5 GeV and 4.6 GeV, respectively. 



obtained as 



M z = (3.91 ± 0.19) GeV. 



(11] 



III. SUMMARY AND CONCLUSION 



By assuming Z c (3900) as a D*D molecular state with quantum numbers I G J P = 1 + 1 + , 
the QCDSR approach has been applied to calculate the mass of the resonance. Our nu- 
merical results are Mz = (3.91 ±0.19) GeV for molecular state. The results are compatible 
with the experimental data of Z c (3900) by BEC and BELLE Collaborations. 
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FIG. 2: The contributions from the pole terms with variation of the Borel parameter M 2 in the 
case of molecular state. The notations a, /3, 7, A and p correspond to the threshold parameters 
/i5= 4.2 GeV, 4.3 GeV, 4.4 GeV, 4.5 GeV and 4.6 GeV, respectively. 
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FIG. 3: The mass of the molecular state as a function of M 2 from sum rule (I10p . The notations 
a, (3, 7, A and p correspond to the threshold parameters t/sq = 4.2 GeV, 4.3 GeV, 4.4 GeV, 
4.5 GeV and 4.6 GeV, respectively. 
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